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ABSTRACT 


The  stability  of  a  viscous  Maclaurin  spheroid  Is  solved 
asymptotically  for  small  kinematic  viscosity,  v  .  It  is  shown 
that,  in  this  limit,  the  frequency  of  oscillation,  n,  with  respect 
to  the  mode  which  becomes  neutrally  stable  in  the  absence  of 
viscosity  at  the  point  of  bifurcation  (where  the  eccentricity, 
e,  of  the  meridional  section  is  approximately  0. 81Z7),  is 

5vnn 

n  =  no  +  i7(i)  +0(v)  (1) 

In  the  foregoing  formula  nQ  denotes  the  frequency  in  the  absence 
of  viscosity,  a  the  radius  of  the  equational  section  and  <p[e)  a 
certain  function  of  e  which  changes  sign  at  e  =  0. 8127  and  is 
positive  for  smaller  values  of  e  .  From  equation  (1)  it  follows 
that  the  Maclaurin  spheroid  is  indeed  unstable  beyond  the  point  of 
bifurcation  when  viscosity  is  present. 
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1.  INTRODUCTION. 

It  is  well  known  that,  for  an  incompressible  fluid  mass  in  a  state  of 
uniform  rotation  and  in  equilibrium  under  its  own  gravitation,  a  spheroid  of 
revolution  is  an  exact  figure  of  equilibrium.  This  is  the  spheroid  of  Maclaurin. 
The  stability  of  these  spheroids  with  respect  to  infinitesimal  perturbations  has 
been  investigated  in  a  number  of  classical  papers  by  Riemann  (’860),  Bryan 
(1888]|  and  Cartan  (1922)  (see  also  Lebovitz,  [1961]).  From  these  investigations 
it  is  known  that,  in  the  absence  of  viscosity,  the  Maclaurin  spheroids  are 
stable  provided  that  the  eccentricity,  e,  of  the  meridional  section  is  less  than 
0.9529;  and  that,  when  the  eccentricity  exceeds  this  limit,  the  Maclaurin 
spheroids  are  unstable  by  oscillations  of  increasing  amplitude  (the  period  of 
oscillation  being,  in  fact,  equal  to  the  angular  velocity  fi  ).  However,  it 
was  known,  even  before  the  upper  limit  e  =  0.9529  was  established  by  Bryan, 
that  the  Maclaurin  spheroid  allows  a  neutral  mode  of  oscillation  at  e  =  0.8127, 
where  the  ellipsoids  of  Jacobi  branch  off. 
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The  occurrence  of  a  neutral  mode  for  the  Maclaurin  spheroids  at  the 
point  of  bifurcation  is  suggestive;  indeed,  it  was  categorically  stated  by 
Thomson  and  Tait  (1883,  p.  333)  that  "the  equilibrium  in  the  revolutlonal 
figure  is  stable  or  unstable  according  as  e  is  less  than  or  greater  than  0.81266.  " 
In  definitely  associating  the  occurrence  of  the  neutral  mode  (for  infinitesimal 
oscillations)  with  the  onset  of  instability,  Thomson  and  Tait  clearly  had  the  role 
of  viscosity  in  mind.  They  state,  for  example,  "if  there  be  any  viscosity, 
however  slight,  in  the  liquid  or  if  there  be  any  imperfectly  elastic  solid,  however 
small,  floating  on  it  or  sunk  within  it,  the  equilibrium  in  any  case  of  energy 
either  minimax  or  a  maximum  cannot  be  secularly  stable:  the  only  secularly 
stable  configurations  are  those  in  which  the  energy  is  the  minimum  with  given 
momentum. 11  All  these  considerations  became,  of  course,  clarified  by  the 
investigations  of  Poincare  which  showed  that,  at  the  point  of  bifurcation  of  the 
Maclaurin  sequence,  it  is  possible  to  deform  the  Maclaurin  spheroid  into  a 
body  of  the  same  angular  momentum  but  of  lower  total  energy.  From  this  last 
fact  one  generally  concludes  that,  should  any  dissipative  mechanism  be 
operative,  then  it  would  indeed  operate  and  induce  the  transition  to  the  state 
of  lower  energy,  i.  e.  instability  will  arise.  In  other  words,  the  belief  and  the 
conjecture  is  that,  were  one  to  investigate  the  stability  of  a  viscous  Maclaurin 
spheroid  (again  for  infinitesimal  perturbations)  then  one  would  find  that  it  becomes 
unstable  (with  an  initial  infinitisemal  amplitude  increasing  exponentially  with 
time)  when  one  surpasses  the  point  of  bifurcation.  This  last  argument  is  usually 
illustrated  by  simple  mechanical  examples  (cf.  Jeans,  [1929],  p.  201).  While 
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the  suggestion  that  the  objects  should  be  unstable  under  the  circumstances 
envisaged  is  very  reasonable,  it  is  not  a  forgone  conclusion  that  motions 
consistent  with  the  equations  of  hydrodynamics  do  exist  which  will  allow 
instability:  examples  are  known  when  a  normal  mode  analysis  proves  stability 
in  circumstances  when  an  energy  argument  would  suggest  instability  (cf. , 
Chandrasekhar,  1961,  S103).  And  finally,  the  energy  argument  provides  no  basis 
for  estimating  the  growth  rates  of  the  unstable  mode  when  viscosity  is  operative. 
For  these  reasons,  we  shall  consider  in  this  paper  the  stability  of  the  Maclaurin 
spheroid  allowing  for  viscosity  in  the  limit  when  it  is  small;  the  assumption 
of  small  viscosity  is  made,  for  the  problem  is  then  tractable  by  a  boundary 
layer  argument. 
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II.  THE  EQUILIBRIUM  STATE 


Let  (x,y,z)  be  cartesian  co-ordinates  in  which  the  z-axis  is  parallel 

to  ,  the  angular  velocity  of  the  spheroid,  and  in  which  the  origin  is  at  the 

center  of  that  body.  The  equation  of  the  surface,  S  ,  of  the  spheroid  is  then 

© 

of  the  form 


~2 

co 

2 

a 


(2.1) 


~  2  2  ~ 

where  <o  =  (x  +  y  )2  .  It  is  convenient,  for  the  later  discussion,  to 
introduce  oblate  spheroidal  co-ordinates  here.  For  this  purpose,  define  k  and 
E  by 

k  =  a2  -  c2,  E  =  c/k  .  (2.2) 

The  following  relations  then  exist  between  a,  c,  k,  E  and  the  eccentricity  of 
the  spheroid,  e  : 

a  =  k(l+  E2)2  ,  c  =  kE  =  a(l  -  e2)2  , 

(2.3) 

k  =  ae,  sin  *  e  =  cot  *  E  . 

Oblate  spheroidal  co-ordinates  (4,  0,  <p)  are  defined  by 

x  =  k(l  4  4  )2  sin  0  cos  <p  , 

2  - 

y  =  k(l  +  4)2  sin  0  sin  f  ,  (2.4) 

z  =  k  4  cos  0  , 
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(0<<«o,  O<0<*,  0  <  <p  <  2  it)  .  They  form  an  orthogonal  curvilinear  system, 
right-handed  in  the  order  (4,  6,  <p),  and  having  scale  factors  h^  hQ  and  h^ 
given  by 


=  k 


1+  42 


hQ  =k(42+nV 


h  *  k(l+  42)^  (1  -  h-2)^ 
<p 


(2.5) 


where  ja  =  cos  0  .  The  surfaces  of  equal  4  are  a  confocal  system  of  spheroids 
and,  in  virtue  of  the  definitions  (2.2),  S  is  the  particular  member 

4  =  E  (2.6) 

of  this  system. 

Since  the  material  of  which  the  spheroid  is  composed  is  supposed  to  be 
incompressible  and  of  uniform  density  p,  the  equation  of  hydrostatic 
equilibrium  is  simply 

£•_  v  -  7  n2  7Z  *  constant  ,  (2.7) 

P  z 

where  p  is  the  pressure,  and  V  is  the  gravitational  potential.  Now  the 
potential  within  the  oblate  spheroid  (2.1)  is  given  by  (see,  far  example,  Routh, 
1922,  Art.  211) 


2 

r"*i-v  "V 


it  Gpa  c 

where  I,  Aj  and  A^  are  constants  given  by 


> 


(2.8) 
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I  ■  /  - 7”^ — 2 - 1  a  ae  sln”1  6  » 

0  (a2  +  u)(c  +  u)2 


00 


du 


^  4  (a2  +  u)2  (c2  +  u)*  (ae)3 


“ — [sin_1e-e(l-e2)2]  ,  (2.9) 


du 


=  J  ~  7  7  *  -  ■»  l  7  i. 

*  0  (a*+u)(c  +u)T  (ae)  (1-e  )2 


2  r  e  ,  -1  . 

—  [ — TI  ~  sin  ®J  • 


Thus,  equation  (2.  7)  may  be  rewritten  in  the  form 

„2 

|  (j  r  * 

•=  constant  -  irGpa  c[(A,  - 


£2  .  /%/  2  2, 

-•  oonctant  —  rr  r*l  /  A  —  — .  —  \  (j  +  A  Z  J  • 

2  3 


(2.10) 


2irGpa  c 

However,  since  p  must  vanish  everywhere  on  the  equilibrium  surface  (2.1)  , 
it  must  be  possible  to  rewrite  equation  (2.10)  in  the  form 


P  _  2  3.  r,  ~2  z2  . 

—  =  vGpa  c  A  [1-  — -  —  ]  . 

r  a  c 


(2.11) 


Comparing  the  expressions  (2.10)  and  (2.11),  we  conclude  that 


2  _ n 

3  - 


)  =  c  A,  , 


(2.12) 


2irGpa  c 


i.  e. ,  by  equations  (2.9)  , 

S22  =  ^p(l-e2)2  [(3-2e2)  sin"1  e-3e(l-e3)*]  . 
e 


(2.12 
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IH.  PERTURBATIONS  ABOUT  EQUILIBRIUM 

Suppose  that  the  equilibrium  solution  of  Section  II  is  infinitesimally 
perturbed.  The  fluid  velocity  u^  in  the  subsequent  motion  will  satisfy  the 
incompressibility  requirement 

div  ^u  =  0  ,  (3.1) 

and  also  the  linearized  Navier-Stokes  equation 

8JL  2 

— ~  +  2fixu  a  -  grad  n  +  vV  u  ,  (3. 2) 

where  v  denotes  the  kinematic  viscosity,  and 

II  =•£-  v  2  .  (3.3) 

P  L 


Thus  u^  must  satisfy  boundary  conditions  on  the  perturbed  surface  S  of  the 
fluid.  These  require  that  the  normal  velocity  of  the  fluid  on  S  should  equal 
the  normal  velocity  of  the  boundary  S  itself.  They  also  require  that  the  normal 
components  of  the  mechanical  stress  tensor  should  vanish  on  S  .  Introducing 
the  rate  of  strain  tensor  R^ ,  defined  by 


8u, 


3x, 


(3.4) 


the  viscous  stress  tensor  is  vpR^,  .  Since  uis  a  quantity  of  the  first  order 

of  smallness,  the  differences  in  the  values  of  R„.  and  R„  on  S  from  their 
J  & 0  Z<P 

value  on  S  are  of  second  order  of  smallness,  and  can  therefore  be  neglected. 

G 
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It  follows  that  the  demands  that  vpR  and  vpR  should  vanish  on  S  require 

±<p 

either  v  =  0,  or  [R^0]g  =  [R^g  =0  >  (3.5) 

©  © 

where  we  have  used  [UJv  to  denote  the  value  of  a  quantity  U  on  a  surface  V  . 
The  remaining  stress  condition  is 

[p-2vPR^]s*0  (3.6) 

and  it  is  again  convenient  to  extrapolate  this  condition  back  to  the  reference 

surface  S  .  For  this  purpose  we  write  p  and  V  in  the  forms 
© 

P  =  Pe+6p,  V  =Ve+  6V  ,  (3.7) 

where  p  and  V  are  the  values  of  p  and  V  in  the  equilibrium  solution 

©  G 

(cf.  eqs.  [2. 8j  and  [2. 11]).  We  also  introduce  ^{x,  t)  the  displacement  of  a 
particle  of  fluid  at  time  t  from  its  equilibrium  position  jx  .  Clearly  we  have 

it  =  8l/dt  »  (3.8) 

to  first  order.  Moreover,  it  is  evident  that 

[pJs  =  [Pe]s^P]S 

e 

=  [pe  +  i*grad  PeJs  +[6pIs 
e  e 


(3.9-) 
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Thus  condition  (3.6)  is  equivalent  to 


[pet|.gradpet6p-2vpRst]s  .  #  >  t3.10, 


or,  since  equations  (2.7)  and  (3.3)  imply  that 


n  =  -  V  -  77  G2  ~  2)  +  (^  -  6V)  =  ^  -  6V  +  Const,  (3.11) 

p  e  L  p  p 


we  have  (absorbing  the  constant  on  the  right-hand  side  of  eq.  [3.11}} 


[®  +  Pn  -  2vPR££ls  =  0  > 


(3.12) 


where 


*  =  p  6V  +  §  •  grad  p 


(3.13) 


The  significance  of  $  may  be  seen  by  considering  the  changes  6W  in 
the  gravitational  energy  W  of  the  configuration 


W 


1_  _  2  r  r  ~  ~ 

=  “  2  Gp  AmV  Ix-x'l  * 


dxdx' 


(3.14) 


consequent  upon  the  deformation  §  .  It  is  evident  from  equation  (3.14)  that 

,  dxdx' 

6V  =  -Gp6  Jv  j 

e  q 


=  -p/v  6Vte)d£  • 


(3.15) 
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On  the  other  hand,  the  work  done  against  the  pressure  forces  in  the  deformation 
^  is,  t©  first  order, 


-/v  6*gradpedx  .  (3.16) 

e  ~ 

Thus  the  sum  of  X 3. 15)  and  (3.16),  namely  /  ®  dx^,  is  the  work  required 

e 

in  order  to  subject  the  equilibrium  configuration  quasi-statically  to  the  £ 
deformation.  In  the  dynamical  problem  studied  herein,  it  will  equal  the  sum  of 
two  parts,  namely,  the  decrease  in  kinetic  energy  and  the  energy  degraded  to 
heat  by  the  action  of  viscosity. 


IV.  GRAVITATIONAL  EFFECTS 


In  calculating  6V,  the  surface  deformation  may  be  divided  into  its  harmonic 
components  Pnm(p)  eim<|P,  and  each  may  be  considered  separately.  We  will 
write,  therefore, 


K  ,  P  eim* 
n’m  n 


(4.1) 


where 


K 


The  gravitational  effect  of  this  displacement  may  be  likened  to  the  addition  of 

a  surface  mass  distribution  to  S  of  density  p[£  ]_  . 

e  * 


The  gravitational 


#37  3 


■11- 


potential  to  which  this  gives  rise  must  be  a  continuous  solution  of  Laplace's 
equation,  and  must  therefore  be  of  the  form 


6V  =  < 


Q  m(i4) 

n  '  _  m.  .  im<p  ..  .  _ 

C.  _  — “ —  P„  (n)  e  ,  4  >  E  , 


n»m  Q  m(iE)  n 


n 


P  m(i4)  . 

n  '  _  m.  .  im <e  ..  .  _ 

c  - - -  P„  M  e  ,  4  <  E  , 


n,mpm  n 
n  '  ‘ 


(4.  3) 


where  Qnm(i4)  denotes  the  associated  Legendre  function  of  the  second  kind 
of  the  orders  n  and  m  indicated,  and  of  imaginary  argument  ,  For  complex 
arguments,  we  adopt  Hobson's  definitions  of  Pnm(z)  (see,  for  example,  Hobson, 
1931,  §  54}: 


o  i  dmP  (zj 
p  m(z)*(Z2-i}*m — sLf 

n  dzm 


m. 


Q  “t*J 


(z2  -  l)*m 


dmC^(z) 


dz 


m 


The  Wronskian  relationship  between  these  functions  takes  the  form 


P 

n 


m,  ,  _  m,  ,  d  _  m,  , 

n  C2>  '  Qn  &  r  Pn  l‘>  * 
z 


(-)m(n+m)l  1 

(n-m)I  X_Z 


(4.4) 


(cf.  Hobson,  loc.  cit. ,  §  148}. 

In  order  to  determine  the  constant  Cn  m  appearing  in  the  solution  (4.  3)  , 
we  apply  Gauss'  theorem  to  an  elementary  volume  containing  the  surface  mass 


distribution: 
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[iT5I5Vl  -[^H6V]E  ' 

%  £=E+0  ?>  ^=E-0  £,=E 


(4.5) 


which,  using  equation  (4.  3),  reduces  to 


C  _ F  jr  Qnm(4)  “  •—  Pnm(iUl  =-4trGpK  .  (4.6) 


LQ  (14) 
n 


Pnm(±U 


By  using  relation  (4.4),  we  finally  obtain 


m. 


r1 • 


(4.7) 


It  follows  that 


4TrGo2i(-nm(n-m)l 
@2  (n+m)l 


nm(iEJ  Qnm(iE)  P„mW  e*"* 


K  P  “(IE)  Q 

n  m  n  v  1  x 


(4.8* 


In  order  to  calculate  [4*]g  ,  we  note  that,  according  to  equation  (2.11), 

e 

grad  pe  =  -  |  GP2  a3  (1-e2)  Agh^  ,  (4. 9) 

where  1^  is  a  unit  vector  in  the  direction  of  increasing  t,  . 

Thus,  by  equation  (4.1)  , 


fegrad  p  ]  =  Gp2  a3(l-e2)A  K  P  m(„) 

e  * 


(4.10) 


Hence 


C*Jc 


.  _  2 

4trGo 

K-l)m(n-m)l 

2 

e 

(n+m)l 

K  P  m(K)irnv 

n  m  n  »*■ 


(4.11) 
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In  the  studies  which  follow,  we  will  be  particularly  concerned  with  the 
case  n  =  m  =  2  .  For  this  harmonic,  we  have 

3 

P,2(iE)  =-3(l+E2)  ,  Q  2 (IE)  =  i^H-E^cot^E-^-^l'*]  (4.12) 

2  2  (1+E2) 

and  $  takes  the  form 

2  x 

[*]s  =  -  Ifif  -  [e(l-e2)I(3+10e2)  -(3+8e2  -  8e4)  sin'1  e]K2  2P2(^e2i*’  . 

e  2e  ’ 

(4.13) 


V.  THE  INVISCID  SOLUTIONS 

Cartan  (1922)  and  others  (see,  for  example,  Bryan,  1888)  have  examined 

the  case  v  =  0  in  great  detail,  and  have  shown  that  all  the  normal  modes  of 

oscillation  can  be  generated  from  polynomial  solutions  of  Poincare's  differential 
* 

equation  : 

[-L2+-V+( l+^)“4]  =0  ,  (5.1) 

9x  8y  s  9z 

where  s  =  9/9t  .  Since  boundary  conditions  (3.5)  are  automatically  satisfied 

when  v  =  0  ,  there  is  only  one  condition  (viz.  eq.  [3. 12])  to  be  applied  on 

2 

Sg  ,  and  this  determines  s  for  the  normal  mode  considered.  Although 

— 

That  Poincare's  equation  arises  may  be  seen  simply  by  solving  eq.  [3.2]  for 
u^  in  terms  of  grad  n  and  substituting  the  result  into  eq.  [3. 1], 
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equation  (5. 1)  is  hyperbolic  when  -4J2  <s  <  0  and  the  boundary  condition 
(3.12)  is  the  Dirichlet  type,  Cartan  was  able  to  establish  that,  even  in  this 
case,  the  polynomial  solutions  are  unique. 

We  can  illustrate  the  solution  simply  for  the  important  case 

2 

n  =  h(x+iy)  ,  (h  =  constant)  ,  (5.2) 


which  evidently  satisfies  equation  (5.1).  From  equation  (3.2)  we  find 


ux  a  _  2h(x±lY)  =  ,2&x±ia  u  =0 

(s-2ii2)  ’  y  (s-2ifi)  ’  z 


(5.3) 


It  should  be  noticed  in  passing  that  this  solution  satisfies 


V  u  =  0  . 


According  to  equation  (5.  3),  u^  is  given  by 


Vs  *  "2  e(1-e2lr  <■ T  V  “lUy  +  "7  uz> 

a  a  c 


(5.4) 


3  2  - 

(s-2iQ) 


h(x+  iy)2 


(5.5) 


Thus,  on  S  we  have 
e 


1  2  2  2  its 

-a  hP  ((a)  e 


2  2  4 

2a  he(l-e V 

3(s-2ifi) 


(H)  e 


2icp 


(5.6) 


(5.7) 
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From  equation  (5.  7),  in  the  notation  of  equation  (4. 1),  we  find 


2  - 

v  _  2a"he(l-e  l2 
2,2  =~  3s(s-2iQ) 


(5.8) 


Thus  drawing  together  equations  (5.  6), ^5. 8),  (4.13)  and  (3.12),  we  find 


2  ..  2i 


+  23s(^-2ii2)  ^^"[e(1-e2)2(3+10e2)-(3+8e2-8e4)sin_1e]  =  0  , 


i.  e. 


_  _  iSfi  / i_q2\I  rwi_Q2\2 


2  „  4.  .  -1 


s(s-2i£2)  =  -  — j*  (1-e  )2  [e(l-e  )2(3  +  10e  )-(3+8e  -8e  )  sin  e]  .  (5.9) 


This  is  Cartan's  result  (loc.  cit.  p.  338).  Writing  it  in  the  form 


s(s-2iQ)  +  Q  s  0  (5.10) 

we  find  that  either  s  =  s  or  s  =  sL  where 

a  b 

2  4 

sg  =  i[«  -  (Q+  £2  )2  ]  =  ina(say)  ,  (5.11a) 

2  - 

sb  =  i[0  + (Q+ n  )  2]  =  ^(say)  .  (5.11b) 

At  the  point  of  bifurcation  of  the  Maclaurin  sequence,  Q  is  zero,  the  root 
(5.11a)  is  therefore  zero,  and  the  spheroid  is  neutrally  stable.  Beyond 
bifurcation,  sq  is  again  non-zero,  and  the  spheroid  is  again  stable.  However, 
if  the  argument  of  Thomson  and  Tait  is  correct  (see  Section  I),  viscosity  should 
introduce  into  root  (5.11a)  a  real  part  which  changes  sign  from  negative  to 

positive  on  passing  through  the  point  of  bifurcation.  This  possibility  is  studied 
in  the  next  section. 


-16- 


#37  3 


Before  leaving  this  solution  it  is  useful  to  compute  the  rate  of  strain 
tensor.  According  to  equations  (3.4)  and  (5.  3)  we  find  that,  in  the  cartesian 
coordinate  system  (x,y,z) 


2h  n  i  °\ 

1  I)  <5-“> 

The  rate  of  strain  is  thus  constant  throughout  the  spheroid.  It  may  be  written 
simply  in  dyadic  form 

%  =-(T^T  a,!,  (5.B) 

where  SL  is  a  vector  whose  (x,  y,  z).  components  are 


&=(1,  i,  0}  , 


(5.14) 


and  whose  (4,  0,  <p)  components  are  therefore 


3,  =  (sin  cos  i)  ei<p  , 


(5.15) 


where 


(5.16) 


By  equations  (5.13)  and  (5.15),  we  find  immediately  that 
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2h  sin2  e2^  2h  slnvli  cos  tii  e2^ 

R4;  =  “  (s-2i£2)  »  R40  =  (s-2i£7) 

(5.17) 

_  2hl  slnjje2^ 

\<p  ~  ~  (s-2i«) 


VI.  THE  EFFECT  OF  A  SMALL  VISCOSITY 

The  solution  obtained  in  Section  V  for  an  inviscid  fluid  can,  in  principle, 
be  modified  in  two  ways  when  the  fluid  has  a  small  kinematic  viscosity  v  . 

First,  the  equations  of  motion  with  v  *  0  may  not  be  satisfied  by  the  inviscid 
solution  so  that  correction  terms  of  order  v  are  required.  In  the  present 
problem,  however,  in  which  we  are  concerned  with  the  solution  (5.  3),  the 
viscous  terms  of  equation  (3.2)  vanish  identically  in  virtue  of  equation  (5.4), 
and  no  such  modifications  of  the  inviscid  solution  are  required.  Second,  the  boundary 
conditions  which  must  be  satisfied  by  the  solution  for  a  viscous  fluid  include 
the  requirements  (cf.  eq.  [3.  5])  that  the  tangential  components  of  R  must 
vanish  at  the  free  surface  S  .  If,  as  in  the  present  problem,  these  conditions 
are  not  satisfied  by  the  inviscid  solution,  then  a  boundary  layer  in  which  viscous 
effects  are  significant  must  be  introduced  in  order  to  make  the  necessary 
adjustments  to  the  tangential  stresses.  We  suppose  that  the  leading  terms  in 
the  boundary  layer  vanish  at  its  inner  edge  so  that  the  inviscid  solution  is  not 
modified  to  the  order  of  these  terms.  However,  in  satisfying  the  conditions  on 
the  tangential  stress  components,  the  condition  on  the  normal  component  is 
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invalidated,  albeit  to  a  higher  order  in  v  .  In  order  to  correct  this,  all  the 
solutions  of  the  inviscid  problem  involving  the  surface  harmonics  for  which 
m  =  2  must  be  excited,  including  the  one  with  n  =  2  .  Equation  (5.9)  is 
consequently  modified. 

2 

For  small  values  of  the  viscosity  (v<<  a  n)  it  is  accordingly  convenient 
to  proceed  as  follows.  Write 

Jio+Ai!  *  ir  =  no  +  n1  ,  *  =  *0  +  *!  >  t6-1) 


where  u„  and  II  „  are  solutions  of  the  inviscid  equations  with  m  =  n  =  2 
but  which  do  not  obey  equation  (5.9),  and  u^  and  ,  are  the  leading  terms 
of  the  modification  due  to  viscosity  which  vanish  in  the  interior  of  the  fluid. 

In  addition  <E>Q  is  the  value  of  $  arising  from  the  perturbation. of  the  free 
surface  due  to  ^u^,  with  a  corresponding  definition  for  .  Then,  from 
equation  (3.12),  we  have 


+  pn  -  2vpr^Js  =  0  > 


(6.2) 


i.  e. 


K * pVs. +  [®i + pni “ 2v|>lWse  1 0  -  l6-3> 


as  the  condition  on  the  normal  stress  to  be  satisfied  by  the  solution.  We  are 

particularly  interested  in  the  harmonic  n  =  m  =  2,  and  so  we  multiply  equation 
2  ™2i 

(6.  3)  by  P2  (|a)  e  and  integrate  over  the  surface  of  the  spheroid, obtaining 
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(6.4) 


where  we  have  used  the  theory  of  Section  V.  Here  y  is  the  contribution  from 
the  second  term  of  equation  (6.3): 


In  the  normal  mode  theory,  it  is  assumed  that  time  enters  the  dependent  variables 
st 

through  the  factor  e  ,  so  that 


^  Ws  =  s  U^S 
6  6 


(6.6) 


Using  equations  (4. 1)  and  4. 13)  and  the  definition  of  Q  (cf.  eqs.  [5.9]  and 
[5.10]),  it  follows  that 


y  32ir 


/  dji  /  d?^  -  t  2sE1Q^h^Ul^  “2vR0cAt 


n  \ 
(l-(x  )  e 


(6.7) 


Our  main  purpose  here  is  to  evaluate  y  . 

An  examination  of  equations  (3.1),  (3.  2}  and  (3. 5)  suggests  that  we  assume 
that  the  boundary  layer  is  of  the  conventional  type  in  which 

u10-O(v*),  u^=0(A  8/8&  =  0(v"*)  , 

8/86=0(1),  8/8?  =0(1)  ,  (6.8) 

so  that 

u1;=o(v),  n1=o(v) 

all  orders  of  magnitude  referring  to  the  dependence  upon_  v  . 


(6.9) 
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The  leading  terms  of  the  8,  <p  and  4  components  of  the  equations  of  momentum 
(3.2)  and  the  complete  equation  of  continuity  (3.1)  may  then  be  written 


s)  Uj0  +  2  0  cos  t|t  u^  =  0  , 


(6.10) 


l~2~2‘s>u ^ eo.*  u,#  -  0  > 
n^  84 


(6.11) 


1  8ni 


(6.12) 


2  8uir  d  2  i  2ik(E2  +  li2)u 

(l+EZ)h  -rf  -  £■  [(1-tx2)2  h  u  ]+ - —k - fL  =  0  . 

4  8&  8>x  8  18  (1+e2)2(1  2 ,2 


(6.13) 


In  all  of  these  equations  h  and  h  are  to  be  evaluated  at  4  =  E  ,  and 

4  o 

are  therefore  functions  of  |j.  alone.  Further,  in  equation  (6. 13),  we  have  made 
use  of  the  property  that  all  dependent  variables  depend  on  <p,  only  though  the 
factor  e  .  Combining  equations  (6.12)  and  (6.13)  we  obtain 

J. 


2  2  2 

=  2  h^  sini)»(l-H.  )uly+-^2  (E  +  k-  )  sin  i|/ (iu^-cos+ uie) 


”  *sE  ^  RlT*2)*  (E2+I»2)^ule]  . 


(6.14) 
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Antlcipating  that  y  is  of  order  v,  we  may  assume  for  the  purposes  of 
simplifying  equation  {6.14}  that  Q=-s(s-2in)  (cf.  eq.  5.10),  whence 

»  2  0+e2)Q 

{(1-H  )[n!  “  2sE  ht  ult^ 

(E2+|jl2)  h  sini{j 

= - f -  [(UjQ-iU^)  (1+COS  iJj)  (s-2li2) 


-{u.  +  iu.  )  (l-cosip)  (s+2iflcosiJj)] 

10  1  cp 

♦  ^e191  “lel  •  l6-15) 

The  values  of  u1Q  and  u^  required  for  the  evaluation  of  the  right-hand 
side  of  equation  (6.15)  follow  from  equations  (6.10)  and  (6.11).  Write 

2  2 

- ~  =  s+2iOcos^  ,  =  s-2ii2  cosijj  ,  (6.16) 

L  “  L  “ 


where  the  real  parts  of  a  and  p  are  both  positive.  The  general  solution  of 

_i 

equations  (6.10yand  (6.11)  which  vanishes  when  v  2  (E-£)  is  large  and  positive 
is  then 


u  _LAe«U-E)  I.B/U-E) 

10  ~  2  A  +  2 


u,  ^iAe^^iBe^  . 
l<p  2  2 


(6.17) 
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The  values  of  the  coefficient  A(n,$>)  and  are  found  from  the  condition 

that  the  tangential  components  of  stress  vanish  at  £  =  E  .  Using  equations 
(5.17),  we  have 


4h. 


■(<*A+(3B) 


2h  slnOi  cosdj  2i<p 
s-2in  8 


(6.18) 


(aA-(3B)  =  - 


21h  sin  Oi 
S-2iI2 


Zi<p 
e  r 


> 


and  consequently 


A  = 


4khE 


2  — 

a(l+ E  )2  (S-2IX2) 


(  oos  4«  — i)  (1— m-  )  e 


2^21* 


(6.19) 


B  = 


4khE 


2  — 

(3(1+ E)2  (s-2ity 


(cosiji  +  1)  (1— jjl  y  e 


It  is  noted  that  the  a  priori  assumptions  (6. 8)  about  the  boundary  layer  have 
led  to  a  reasonable  solution  which  may  easily  be  shown  to  be  consistent  provided 
neither  a  nor  p  vanishes,  i.  e.  if  we  write 

s  =  2ifi  cos  ,  (0<  4<0  <  II)  *  (6.  20) 


provided  we  exclude  the  neighborhood  of  the  singularities  ij/  =  ij<  and  vj>  =  n  -i|iQ 

* 

vj>g  is  real  for  the  range  0.  4285  <  e  <  0.  9529  ...  of  eccentricity.  Since  this 
range  contains  the  point  of  bifurcation,  the  singularities  defined  by  equations 
(  6.  20)  are  germane  to  our  discussion. 
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These  zones  require  special  treatment  which  we  shall  postpone  until  the  next 

L 

section.  We  merely  observe  here  that  their  contribution  to  y  is  of  order  v*  . 

Substituting  the  results  (6.17}  and  (6.19)  into  the  first  term  appearing  on 
the  right-hand  side  of  equation  (6.15)^  integrating  with  respect  to  t,  from 
-oo  to  0,  and  noting  that  n,  and  vanish  in  the  interior  of  the  fluid,  we 
find  that 


e 


2hL  hv  (1+E2) *  (1-p2)*  sin  ijj  eZi<p 
m  kE(s-2ifi) 


[(1+  COS  ij l)Z  +  (1  -  COS  »p)2] 


hf  s-2ifl)  8 
2E  8|ji 


[(l-n2) V  +  nV  ule]dS  . 


(6.21) 


The  determination  of  y  is  completed  by  using  the  form  of  R 
equation  (5.17) 


given  by 


r2  R  i  -  .  41*vE2  (1-/)  e21’’ 

1  °“V  (EVWiO)  ’ 


(6.22) 


and  adding  it  to  equation  (6.  21)  giving 

p(i-h  )[n2  -  2sE^  "2vRo;4^S£ 

0 


8hvp  2  2i<p  t  okfs-2inl  8  n  2,f-  2  2 A  , 

“(s-2in)^  +  j  w  2E  8|a  K1”*1  ^  (e  +1*  )  u10]  . 


(6.23) 
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Carrying  out  the  integration  (6.  5),  assuming  that  the  orders  of  integration 
with  respect  to  £  and  pi  can  be  inverted,  we  obtain 


lOhvo 

Y  =  3(s-2in) 


(6.24) 


and  hence,  from  equation  (6.4), 


“L 


S(S-2W)  a2(s_ 


=  0 


2  in) 


(6.25) 


is  the  modified  dispersion  relationship  for  the  harmonic  n  =  m  =  2,  neglecting 
powers  of  v  above  the  first.  Let  one  of  its  roots  be 


s  = 


VS1 


(6.26) 


where  sQ  is  one  of  the  roots  (5. 11)  and  s1  is  of  order  v  .  On  substituting 
into  equation  (6. 25)  and  neglecting  all  powers  of  s^  above  the  first,  we  find 
that 


5v8. 


5vs, 


Sl  =  - 


a  (s  -2i«) 


a2Q 


(6.27) 


2 

For  small  values  of  fi,  equation  (5.11)  implies  that  Q  =  -  s^  .  Thus, 

2 

according  to  equation  (6. 27)  s^  =  -  5  v/a  ,  in  agreement  with  the  result  of 

Lamb  (1881)  for  the  oscillations  of  a  slightly  viscous  sphere  (see  also 

2 

Chandrasekhar,  1959).  Since  sQ  is  always  negative  up  to  the  point  of 
ordinary  instability  (e  =  0.9529),  the  stability  of  the  slightly  viscous 
spheroid  depends  on  the  sign  of  Q  .  Now  Q  is  positive  between  e  =  0  (the 
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sphere}  and  e  =  0. 8127  (where  the  Maclaurin  sequence  bifurcates)  and  is 
negative  between  this  value  of  e  and  e  =  0.9529.  It  follows  that,  beyond 
bifurcation,  the  value  of  s^  is  real  and  positive  so  that  oscillations  of  the 
Maclaurin  spheroid  in  which  the  harmonic  n  =  m  =  2  are  excited  increase 
exponentially  in  amplitude  with  time  in  agreement  with  the  conclusion  based  on 
the  energy  argument  of  Section  I.  Presumably  once  this  instability  sets  in, 
the  configuration  moves  towards  the  lower  energy  Jacobi  state  with  the  same 
angular  momentum. 


VII.  THE  SINGULAR  CIRCLES 

The  theory  of  the  boundary  layer  discussed  in  the  previous  section  failed 
in  the  neighborhood  of  the  circles  of  latitude  +  =  ^0,  and  iJj  =  II-iJ>0,  where 
4<0  is  defined  by  equation  (6. 20).  One  still  expects  the  boundary  layer  to  be 
vanishingly  thin  and  the  velocities  in  it  to  be  vanishingly  small  as  v  -►  0,  but 
the  hypothesis  (6.  8}  on  which  the  theory  is  based  needs  modification.  Let  us 
concentrate  attention  on  the  zone  R+  in  the  neighborhood  of  4'  =  and 
suppose  that,  instead  of  assumption  (6.  8), 

=  0(v_P),  9/8^  =  0(v"q)  ,  (7.1) 

where  1/2  >  p  >  q  >  0;  thus  we  expect  the  boundary  layer  to  be  thicker  in  R+ 
and  rapid  changes  to  occur  in  the  ^-direction  hut  that  the  most  rapid  changes 
still  occur  in  the  £ -direction.  It  will  be  shown  that  p  =  2/5,  q  =  1/5  leads 
to  a  consistent  picture. 
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Further  we  shall  assume  that 


uie  =  °(vr),  %  =  0(vr)  , 


(7.2) 


where  r  >  0  .  The  equation  of  continuity  (6.13)  may  now  be  written  as 


6u,„  ,,  2, j  8u 


-L  _  Chi h  i 

94  he 


16 


dp. 


=  0 


(7.3) 


where  h^  and  h0  are  now  constants,  so  that  q  \  .  The 

equation  of  momentum  in  the  6-direction  reduces  to 


V  aZui8 

— r  - —  =  su.a  -2  Si  cos+  u 

h  2  a;2  10 


(l-u2l*  8ni  . 

1 V  h  e  dp. 


(7.4) 


notice  that  in  R+,  su1Q  =  2 O  cos  ij;  u^  by  equations  (6. 16)  and  (6.17) .  It 
follows  from  equation  (7.4)  that  =  0(v^  +  r  2p+q)  .  The  equation  of  momentum 
in  the  <p-  direction  reduces  to 


a2ui 

77  ~T£=  s%  +  2  0(008+  u10+simliu1;)  , 
h^  64 


(7.5) 


the  pressure  terms  being  of  smaller  order  than  the  other  terms  of  equation  (7.  5) 
since  dTl^/d<p  =  2iII^  .  The  equation  of  momentum  in  the  4-direction  is 


1  0II1 


^sr=2nvin*  > 


(7.6) 
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the  viscous  terms  being  negligible  because  an^/8^  »  811^  /8^  in  virtue  of 

assumption  (7.1),  and  u,  <<  u.  in  virtue  of  assumption  (7. 2)  and  equation 

lg  19 

(7.3).  It  is  to  be  particularly  noted  that  this  equation  is  identical  with 
equation  (6.12),  and  that  both  the  equations  used  to  derive  the  basic 
result  (6.14)  hold  also  in  R  +  .  Consequently,  our  aim  here  is  to  find  the 
contribution  of  R+  to  equation  (6. 14).  If  we  now  eliminate  and  u^  from 
equations  (7. 3)  to  (7.6),  and  write 


V 


10 


-  iu 


lf> 


(7.7) 


we  find  that 


v 


8Z+  hr  2  i 

(s-21X2cosi^)  T7~  ~  2ii2  sinip  r~(l-n  )2 
ol  n„  t 


(7.8) 


Similarly,  if 


Z-  =  U10  * 


83Z 


94' 


8Z 


(s+2in  cosi}>)  -TT-  +  2itt  sinifi  ^(l-(i2)2 


!  8Z 


(7.9) 


(7.10) 


Thus,  in  equation  (7.  8),  the  last  term  on  the  right-hand  side  may  be  neglected 
except  in  R+,  and,  in  equation  (7.9),  the  last  term  may  be  neglected  except 
in  R_,  the  neighborhood  of  9  =  it  -  .  Concentrating  attention  on  R+  , 

let  us  write  9  =  »}j  +  ¥  where  *  is  small.  Then  equation  (7.8)  simplifies 
to 
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v 


2iO  sin  4<q 


210  . 
kE  \ 


sin 


(7.10) 


From  a  comparison  of  the  orders  of  magnitude  of  the  three  terms  composing 
equation  (7.10),  and  by  using  the  assumed  properties  (7.1),  we  have 

a  =|  ,  b=^  .  (7.11) 


Outside  R+ 

Z+  =  BeP^"E* 


(7.12) 


by  equations  (6.17).  Consequently  on  the  edge  of  R+  where  is 

small  but  v  *  is  large, 


V 


X 

g  v2 

7“T  exp[{2iOh  2  sirnj/  f}2(C~E)]'  , 

l*i2  4  0 


(7.13) 


where  is  a  factor  of  order  unity.  Comparing  the  orders  of  magnitude,  it 
follows  that 


Z+  =  0(v2A)  (7.14) 

in  R  .  Moreover  9Z  /9  =  0(1)  in  R  ,  as  it  must  be  since  it  is  required  to 

T  T  ^  T 

play  a  part  in  the  cancellation  of  the  tangential  stresses.  Although  it  has  not 
been  explicitly  evaluated,  the  determination  of  Z,  in  R,  seems  to  be  a 

T  + 

properly  posed  problem  since,  for  not  only  do  we  have  the  behaviors  (7.13)  to 
determine  the  relevant  solution  of  equation  (7.10)  at  the  edges  of  R+,  but  we 
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also  have  Z+**0  as  £  -*•  and,  finally,  8Z+/8^  =  known  stress  on  &=E  . 

Referring  back  to  Section  VI,  the  presence  of  the  singular  lines  led  to  two 
apparent  difficulties  peculiar  to  the  boundary  layers  studied  here.  First,  it 
was  not  clear  that  the  orders  of  integration  and  differentiation  can  be  inverted 
as  was  done  in  evaluating  the  integral  (6.  5}  from  the  result  (6.23).  Now  that 
the  behavior  of  Z,  in  R,  and,  by  implication,  of  Z  in  R  is  seen  to  be 
satisfactory,  we  can  justify  the  inversion  by  dividing  up  the  interval  -l<.p<.l 
into  three  regions  by  separating  out  R_  and  R+  .  In  none  of  these  three 
regions  does  the  inversion  present  any  difficulty.  Second,  a  question  arose 
concerning  the  contribution  from  R+  and  R_  to  y  via  the  first  term  on  the 
right-hand  side  of  equation  (6.15).  The  contribution  from  R+  is  easily  seen 
to  be  of  the  same  order  in  v  as 

///  d?d£d¥  Z+(s-2i ft  cos  +) 

i*  6*  dS 

///d?d;d*z+*  (7.15) 

and,  using  the  results  (7.11)  and  (7.14),  this  integral  is  of  order  v*/*  .  Thus 
we  see  that,  although  the  velocities  in  the  neighborhood  of  R+  are  larger  by 
a  factor  of  v  V»o  than  those  elsewhere  in  the  boundary  layer,  their  contribution 
to  y  is  smaller  by  a  factor  of  v  V*  .  A  similar  remark  applies  to  R  . 
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